Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



o 

Q 

o 



> 
OS 

(N 

o> 

OS 

0,; 

6 

C3 



13 



JET ACCELERATION BY TANGLED MAGNETIC FIELDS 
Sebastian Heinz 1 

AND 

Mitchell C. Begelman 2 ' 3 

JILA, University of Colorado and National Institute for Standards and Technology, Boulder, Colorado 

80309-0440 4 
Draft version February 1, 2008 

ABSTRACT 

We explore the possibility that extragalactic radio jets might be accelerated by highly disorganized 
magnetic fields that are strong enough to dominate the dynamics until the terminal Lorentz factor is 
reached. Following the twin-exhaust model by Blandford & Rees (1974), the collimation under this 
scenario is provided by the stratified thermal pressure from an external medium. The acceleration 
efficiency then depends on the pressure gradient of that medium. In order for this mechanism to work 
there must be continuous tangling of the magnetic field, changing the magnetic equation of state away 
from pure flux freezing (otherwise conversion of Poynting flux to kinetic energy flux is suppressed) . This 
is a complementary approach to models in which the plasma is accelerated by large scale ordered fields. 
We include a simple prescription for magnetic dissipation, which leads to tradeoffs among conversion of 
magnetic energy into bulk kinetic energy, random particle energy, and radiation. We present analytic 
dynamical solutions of such jets, assess the effects of radiation drag, and comment on observational issues, 
such as the predicted polarization and synchrotron brightness. Finally, we try to make the connection 
to observed radio galaxies and 7-ray bursts. 



1. INTRODUCTION 

Many extragalactic radio jets move with bulk Lorentz 
factors r > 10, as evidenced by very short variability 
timescales, superluminal proper motion of jet features, and 
dilated particle lifetimes, yet the process that actually ac- 
celerates the material is not well known, and neither is the 
mechanism that collimates the outflow. 

One of the first serious models for the large scale dynam- 
ics of extragalactic radio jets is the 'twin exhaust' model 
(Blandford & Rees 74, BR74 throughout the rest of the 
paper). In this model, relativistic particle pressure pro- 
vides the bulk acceleration via conversion of internal to 
kinetic energy. The collimation comes from confinement 
by an external medium, the pressure of which is very likely 
stratified in the gravitational field of the central black hole 
and the host galaxy and cluster (either in a hydrostatic 
equilibrium configuration or as a wind). Pressure balance 
then dictates the evolution of the jet bulk Lorentz factor 
r, and the jet diameter R. However, cooling processes for 
particles with highly relativistic random motions (neces- 
sary to produce outflows with large bulk Lorentz factors 
r) are very efficient, thus competing with and probably 
disabling bulk acceleration. 

The problem of radiative losses is similarly limiting in 
the case of radiatively accelerated jets (Phinney 1982). 
Furthermore, these jets rely on the presence of a strong 
point-like radiation source (the terminal Lorentz factor 
is limited by the solid angle subtended by the radiation 
source due to relativistic aberration). The same Inverse 
Compton (IC) scattering effects that lead to the bulk ac- 
celeration lead to even stronger radiative losses. In fact, 



radiation drag can actually decelerate the jet if the exter- 
nal radiation field is isotropic, which leads to both random 
energy losses and kinetic energy losses by the same process. 

Radiative losses are much less limiting if the energy is 
stored in the magnetic field. Organized magnetic fields 
have been suggested to provide both acceleration and col- 
limation (e.g., Blandford & Payne 1982, Li, Chiueh, & 
Begelman 1992), however, such solutions, in which the col- 
limation is provided by toroidal field, seem to be hampered 
by instabilities (Begelman 1998). Furthermore, it is diffi- 
cult to achieve collimation on the basis of magnetic ten- 
sion alone without contribution from an external pressure 
source (Begelman 1995, B95 hereafter). 

Polarization measurements show that the magnetic field 
in jets is probably not well organized — the polarization 
is generally well below the maximal value of ~ 70%; only 
in knots does the polarization tend towards this value 
(note, however, that interpretations of polarization mea- 
surements are often ambiguous, since different field ge- 
ometries can sometimes lead to the same net polarization) . 
This argues for the presence of largely unorganized, chaotic 
fields, which could easily account for the high polarization 
measured in the knots if they are interpreted as shocks, 
compressing the field in the shock plane (Laing 1980, B95). 
This goes hand in hand with the fact that the field pro- 
duced in the disk by dynamo processes is expected to be 
highly chaotic. Since the conditions in the jet will likely 
be controlled by disk physics, we should expect the same 
statement to be true for the magnetic field in the jet at 
least close to the disk. These arguments led us to investi- 
gate the dynamics of jets containing large amounts of such 
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disorganized magnetic fields. 

The rate of acceleration in a jet propelled by internal 
(isotropic) particle pressure in an external pressure gradi- 
ent is limited to T a p C xt -1 ^ 4 (BR74), which means that 
the acceleration to bulk Lorentz factors of T ~ 10 — 100 
would occur over length scales > (1000) r g for external 
pressure gradients p oc z~ 2 . One might think that an 
anisotropic pressure in the form of chaotic magnetic fields 
could increase the rate at which the jet is accelerated, if 
the excess momentum flux is oriented along the direction 
of the jet. We will show that under a given set of simple 
assumptions the rate of acceleration is actually the same as 
in the classic case considered by BR74, i.e., T cx Pcxt -1 ^ 4 - 

It is unlikely that the magnetic field evolves without 
some form of dissipation, especially if it is highly unorga- 
nized (reconnection is a diffusive process, so strong gradi- 
ents in the field, as are present if the field is highly tangled 
on small scales, will likely lead to increased dissipation). 
These loss processes can compete with the efficiency of 
bulk acceleration by removing energy from the flow reser- 
voir. We will investigate the effects such a tradeoff might 
have on the dynamics and appearance of jets. 

We will lay out the simplifying assumptions going into 
this model in §|^. Section y contains a brief discussion 
of allowed solutions and the sonic transition of such so- 
lutions. In §0 we will discuss self-similar and asymptotic 
solutions ana present some full analytic solutions, includ- 
ing a simple estimate of the effects of radiation drag. Sec- 
tion [5] contains a discussion of the aforementioned tradeoff 
between dissipation and acceleration, some predicted ob- 
servational consequences, and an outlook on applications 
of this model to radio galaxies and 7-ray bursts. Finally, 
we will summarize our results in 

2. THE MODEL 

The model we are employing here is closely related to 
(and an extension of) the 'twin-exhaust' model put for- 
ward by BR74. We adopt a similar scenario under which 
the jet is launched into a stratified external medium. In 
the case we are considering, turbulent and highly disor- 
dered magnetic fields dominate the internal energetics of 
the jet. 

We have tried to illustrate the overall picture we are em- 
ploying in Fig. [l| Interstellar magnetic field is advected 
inward by the accretion disk. Turbulent shear then am- 
plifies the field and tangles it up (dynamo action). This 
effect will grow stronger with decreasing distance to the 
black hole. Eventually, regions of very high field strength 
will develop. Due to their buoyancy they will accelerate 
away from the black hole, forming an initial outflow. This 
outflow is then collimated by the pressure of the external 
medium. The jet channel is constrained by pressure bal- 
ance, i.e., the jet will expand or contract in such a way 
that an equilibrium solution is set up for which the flow 
is stationary. As the flow expands, we assume that mi- 
cro instabilities and turbulence constantly rearrange the 
field. As in the pure particle pressure case, the flow can 
go through a critical point, where the radius R has a mini- 
mum and beyond which the flow will become self-similar (if 
the external pressure itself behaves self-similarly with dis- 
tance to the black hole) before the rest mass energy starts 
dominating the inertia, at which point the jet will reach a 



terminal Lorentz factor Too . Along the way, the field might 
dissipate energy via reconnection-like processes, and radi- 
ation drag might alter the dynamics. 

We assume that no energy or particles are exchanged 
between jet and environment, except for radiative losses. 
However, the momentum discharge (the jet thrust Q) need 
not be conserved along the jet. By assuming that quanti- 
ties like p', U'i = (B\ )/8tt (where a prime denotes that 
the quantity is measured in the comoving frame) do not 
vary significantly across the jet we simplify the analysis 
to a quasi- ID solution. We ignore effects of shear at the 
jet boundaries. We assume throughout most of the paper 
that the advected matter is cold (i.e., enthalpy density 
b! w n'm par ticieC 2 ). We are looking for stationary flow 
along the jet (i.e., far from the terminal shock), enabling 
us to drop time derivatives. Finally, to make this quasi ID 
treatment possible, we will need to make the assumption 
that the jet is narrow, which in our case implies that the 
opening angle is small compared to the beaming angle, 
i.e., dR/dz <C 1/r. As we will later show, this also im- 
plies that the jet is in causal contact with its environment 
(as required by the assumption that the jet is in pressure 
equilibrium with the surrounding medium). 

2.1. Treatment of Magnetic Field 

We use cylindrical coordinates (r, <f>, z) with the z-axis 
oriented along the jet axis. The flow velocity is not aligned 
with the z-axis for r > (the jet expands). We assume 
that the sideways velocity is small compared to v z but non- 
vanishing, i.e., the flow is well collimated. The magnetic 
field is expressed in a different basis, since the standard 
basis vectors e r and e z are not orthogonal in the comoving 
frame for r ^ 0. One axis of this new basis is aligned with 
the local velocity vector, e||. The second unit vector 
is coincident with the 0-unit vector of the lab coordinate 
system. The third unit vector e ro is obtained from the 
cross product of the other two. In the comoving frame we 
have 

B' =B' ZB e a + B'^ + B' l] e ]] . (1) 

As mentioned before, the ID approximation is only pos- 
sible if the opening angle is small compared to the beam- 
ing angle. This is because the Lorentz factor will not be 
nearly uniform across the jet otherwise. The assumption 
that dR/dz <C 1/T simplifies the equations of motion sig- 
nificantly. 

Following B95, who investigated similar jets in the non- 
relativistic limit, we assume that the magnetic field is 
highly disorganized. In the comoving frame, averages over 
the individual components and cross terms vanish while 
the energy density in the individual components is not 
zero: 

(B\B'j)=0, fori^j, (B\) =0, 

(B\ 2 ) = 8nU! + 0. (2) 

Lorentz transformation of the field to the lab frame (and 
to the cylindrical coordinate system aligned with the jet 
axis) yields 

B = (^rS' ro + e r + TB'^ 

+ f-S'ii - -rsO e z . (3) 
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Some of the components are now correlated. The electric 
field in the lab frame is 



E 



v.TB'^. (4) 



2.1.1. 



Magnetic Equation of State 

Without the presence of turbulent rearrangement of the 
field, flux freezing would govern the behavior of the indi- 
vidual components. If we assume the presence of turbulent 
mixing between the different field components, we might 
expect the field to follow a modified evolution according 
to 



dB'i 



dB's 



j 



dTv 



d(Tv) 



dR 



dR, (5) 



where the subscript ff denotes the value the derivative 
would take under flux freezing, and ay and are arbi- 
trary mixing coefficients. Based on this picture we there- 
fore choose the following convenient ad-hoc parametriza- 
tion of the field evolution with Lorcntz factor T and jet 
radius i?, including rearrangement: 



B'm (X B'a, 



(vry 2+ ^ r- 2+ ^, 



5' 1, 2 oc {vTf 3 R- 4+ ^ 



(6) 



This is the magnetic equation of state we use. In the case 
of pure flux freezing, //; = for all i. In the case of a 
completely isotropic field we have [i\ = 2/3, \X2 = —2/3, 
/i3 = —4/3, [i = 4/3. Note that this prescription is 
still fully general [until we make some limiting assump- 
tions about the /Ui(r, z)\. Since the re-arrangement process 
mixes the perpendicular and parallel components of the 
field, we would expect that the field behavior is changed 
from flux freezing in such a way that the coefficients /ii 
are bracketed by the values they take in the case of flux 
freezing. Since the case of a purely isotropic field must 
be included in our analysis, it is clear that this condition 
requires that fi2 < < \i\ and pt^ < < fi2- 

We define two quantities to characterize the anisotropy 
of the magnetic pressure: 



u'w-u'± 

U'\\ + U'± ' 



5 = 



UU + U' 



(7) 



where U'± = U' m + U'<p. Thus, the magnetic field is 
purely perpendicular for £ = — 1, and purely parallel for 
f = 1. The perpendicular component is purely radial for 
5 = — 1 and purely toroidal for 6 = 1. The field is per- 
fectly isotropic for ( — —1/3 and 5 — 0. It is obvious 
from equation (^|) that 8 is constant for any combination 
of parameters, since U' m oc U' <f> by assumption. 

While this paramctrization alone is rather unrestrictive, 
we can limit it to a one parameter family by assuming that 
the /ii are constants under any possible variation of T and 
R, and that the rearrangement process does not change 
the total comoving energy density in the magnetic field. 
(Otherwise the same process would have to act as an en- 
ergy sink, since we assume that the magnetic field is the 
dominant term in the internal energy budget. We would 
therefore be de aling with a dissipative process, which we 



will address in §2.1.2.) We can solve for /ij in terms of £ by 



density U' = U\ behave the same as it would following 
flux freezing: 



dU' = U'±_ [(mi - 2) d(Tv) + (/x 2 - 2) dR] 
+U' {1 [» 3 d(Tv) + {Li 4 -4)dR] 
= -2U'± (d(Tv) + dR) - AU\\dR 



(8) 



for arbitrary d(Tv) and dR. Constancy of any of the /ii 
then implies constancy of C an d substitution of £ from 
equation (^) yields 



Ui cx (Tv)^ 1 R- 3 ^ 

Mi = l + C, M2 = -l-C) 
M3 = C-1, (H = l-C 



(9) 



which includes the isotropic case, where the magnetic field 
behaves like a relativistic gas, for which ( = —1/3, (5 = 0. 

It turns out that one can find special analytic solutions 
with constant that satisfy equation (§]) without the re- 
quirement that C, be constant (see § ^) . For these solutions 
the rearrangement process conserves the comoving mag- 
netic energy density only under the variations in T and R 
allowed by the Bernoulli equation [i.e., d(Tv) and dR in 
equation (g) are not arbitrary] . The only condition on the 
parameters fii for such a solution is that ^1/^2 — M3/a*4- 
These solutions arc limited to the self-similar range, where 
the jet is dominated by magnetic pressure. Once they ap- 
proach the terminal phase (i.e., p' J> U'), the parameters 
Hi must vary with z. For the rest of the paper we will as- 
sume that equation (^) holds unless indicated otherwise. 

2.1.2. Dissipation of Magnetic Energy 

It is unlikely that the tangled magnetic field evolves 
without any dissipation of its energy (e.g., via reconnec- 
tion). We thus include a simple, ad hoc prescription of 
magnetic energy losses. We base our parametrization on 
the idea that the magnetic field is always in a nearly force- 
free equilibrium. However, it is impossible to maintain 
perfect force-free conditions everywhere and as the jet ex- 
pands in either direction, the field responds by rearrange- 
ment between the different components (eq. and by 
dissipation of some of its energy. We therefore assume 
that the dissipation rate is roughly proportional to the 
divergence of the velocity in the comoving frame: 



dU'i 

dT / diss 

or, in the lab frame 



-AC/MV-v'l 



dU'j 

dz 



-AU'i 



d: 



In (TvR 2 



(10) 



(11) 



fixing cither T or R and demanding that the total energy 



This ansatz can easily be generalized to different Aj for 
different field components (e.g., if the Alfven velocity fac- 
tors into A). For a more realistic dissipation model see the 
Appendix. 

We will assume that the dissipated energy goes into 
isotropic particle pressure, which is then either (a) radi- 
ated away immediately as isotropic radiation in the comov- 
ing frame or (b) accumulated until the jet reaches a state 
of equipartition between particle pressure and magnetic 
field. 
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2.2. Equations of Motion 
We write the relativistic continuity equation as 

p'v z TR 2 = const. 



(12) 



The energy and momentum equations are given by 
T al3 .p = 0, (T is the stress-energy tensor, separable into 
a matter and an electromagnetic part). In the absence 
of gravity, this reduces to T a ^ ^ — 0, which will be suf- 
ficient for the analysis through most of this paper since 
most of the acceleration will likely take place at distance 
z ^S> r g = GM/c 2 , where M is the mass of the central 
black hole. It turns out, however, that gravity is impor- 
tant in discussing the critical points of the jet, in which 
case we approximate the covariant derivative by a Newto- 
nian potential (-r° 30 = r 3 33 = r 3 00 w r g /z 2 ). We will 



comment on the accuracy of this approximation in §3.1 



Since there is no energy exchange between the jet and 
the environment, and using the expression for the electro- 
magnetic field measured in the lab frame from §2.1, we 



can write the energy equation as T 0a _ a — (neglecting 
gravity). We integrate the equation over a cross-sectional 
volume of the jet and convert it to a surface integral using 
Gauss's law. The contribution from the sidewall is zero, 
giving 

„'„2 



r 2 (p'c 2 + V) v z ttR 2 



4tt 

L = const 



(13) 



(where the B'\ are now averaged quantities). Divid- 
ing equation jl^ ) into equation ([b|) gives the relativis- 
tic Bernoulli equation. In the more general case including 
radiative losses and gravity we have 



+2^T 2 v z R 2 (p'c 2 + Ap' + 2U' X ) + Srad = 0, 



(14) 



where S ia d is the energy lost to radiation leaving the jet. 
We have to make some assumption about the form of S ra d, 
i.e., th e amo unt of energy radiated away [cases (a) and (b) 
from §E.1.2j. 

The z-momentum flux Q can be calculated in much 
the same way (integrating T 33 across a jet cross-section). 
Since the jet can exchange z-momentum with the environ- 
ment, the momentum discharge need not be conserved, 
however. Dropping terms of order v r 2 , the integration 
yields 



Q= dT 

J A 



33 



T 2 v 2 ttR 2 (p + 4p'/c 2 ) + irR 2 p' 
ttR 2 [r 2 (1 + v 2 ) U' ± - U\] . (15) 



The sideways momentum equation is given by T la :a — 
VrT 0a - a = 0. The condition that the solution be stationary 
(i.e., dR(z)/dt = 0) gives the pressure balance condition 
between the jet and its environment. We assume that the 
internal structure of the field adjusts to maintain the given 
cross-section. Since we assume that v T <C v z , the internal 



variation will be sufficiently small, ~ (it) Jz' ^° J us ~ 
tify the assumption of uniformity (note: this assumption 
is only satisfied if the jet is in causal contact). We are thus 
only interested in the pressure balance condition at the jet 
walls, r = R, which gives 



Pe xt = U^ + U' {l -U^+p'. 



(16) 



Note that = directly at the jet boundaries, since 
the magnetic field is assumed not to penetrate the contact 
discontinuity. However, since interior pressure balance de- 
mands that U'q + U!i — U^j + p' be constant, we can set 
U'a — U^, — const, and substitute it for U'i at r = R, which 
gives equation (jig). 




Fig. 1. — Cartoon of the general picture employed in this article. 
Tangled magnetic field is generated in the disk, advected inward by 
the disk flow, and accelerated away from the black hole to form an 
initial outflow. Under suitable conditions, this outflow is then colli- 
mated and accelerates away from the core, keeping pressure balance 
with the thermal pressure provided by the jet environment. 



3. DYNAMICAL SOLUTIONS 

Before we start analyzing the equations presented above, 
it is worth noting that in the case of a cold (p' = 0) jet and 
a magnetic field following pure flux freezing (pi = fj,2 = 0) 
the only possible solution to equation ( |l3| ) far away from 
the core (i.e., z 3> r s ) is T(z) = const., i.e., the jet expands 
sideways to satisfy pressure balance, without accelerating. 
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This is because both the kinetic energy flux and the Poynt- 
ing flux do not vary with R, while they do vary with v, so 
that equation (|l3| ) becomes an equation of T only. Thus, 
fixing the total jet power L fixes T. While a scenario like 
this might explain the coasting phase of the jet (where no 
more acceleration occurs), it cannot account for the initial 
bulk acceleration we are looking for. 

Note that this is different than the case of anisotropic, 
relativistic particle pressure in the absence of isotropiza- 
tion (i.e., simply under adiabatic behavior of the individ- 
ual components). In that case, the components scale like 
p z cx (Tv)~ 2 R~ 2 , p r ocp0 oc (Tv)^ 1 R~ 3 . We might expect 
a behavior like this for a relativistic turbulent pressure 
term. The sideways pressure is simply p r . At relativistic 
speeds, the solution approaches the one found by BR74, 
r cx R cx p -1 / 4 . Thus, unlike in the magnetic case de- 
scribed in this paper, it is generally possible to accelerate a 
jet with anisotropic particle pressure without making any 
arbitrary assumptions about the randomization process. 
This is simply because only the perpendicular component 
of the field, U' ±, contributes to the Poynting flux, while all 
components of the pressure enter equation ([l3]) , which in- 
troduces a dependence on i?, making a solution T =/= const. 
possible. For a magnetically dominated solution to exist, 
on the other hand, we need a field rearra ngement process 
at work, such as was described in §2.1.1. But even under 
such favorable conditions, a proper, accelerating solution 
is not always guaranteed. 



3.1. Critical Points 

Since the jet will likely be injected with sub-relativistic 
speed, the question arises as to where the jet crosses possi- 
ble critical points and at what velocity. If the jet is injected 
at large distances from the central black hole, we can ne- 
glect gravity; if it is injected close to the hole we will have 
to include at least a phenomenological gravity term. 

As a first step, we will set M = artificially (still as- 
suming the presence of an external pressure gradient) and 
neglect dissipation (S ra d = 0). Equation (|lj) gives 



r 2 [(p'c 2 + 4p' + 2U' ± ) v 2 
+ 4(l- 7ad )p' + 2C[/'J 



dv 

vdz 



dR 



+ {4(2- 7ad )p' -2(l + QU'±} — = 0, (17) 

Rdz 

where 7 a( j is the adiabatic index of the particles, 7 a d = 
dhip' /dhip' . This equation has a critical point when the 
expression in square brackets vanishes. At such a point, 
the jet cross-section must satisfy dR/dz = 0, i.e., the jet 
must go through a nozzle, the position and cross section 
of which are determined by the dynamics of the flow. Fol- 
lowing the notation of BR74, the velocity at which that 
happens is 



(7ad - 1) - 2U' ±i ,C 
p'+c 2 + 4p'+ + 2U>^ ' 



(18) 



where the subscript * indicates that the quantity is eval- 
uated at the critical point z*. Since for a magnetically 
dominated jet lim^o- c* = 0, the critical point exists only 
for £ < 0. For a purely isotropic field, where £ = —1/3, 



c* reduces to the sound speed of a relativistic gas with 
7ad = 4/3, = VT73. 

Locally we can always write p cx t cx z~^, thus we define 

£ = -dlnpext/dlnz. (19) 



We can then substitute the pressure balance condition (16) 
into equation ([l7]) in the limit p'c 2 + Ap' <C U' and elimi- 
nate R, which yields 

r 2 [(3 + c)« 2 + (i + 30] 7 ^-(i + C)|. (20) 

This equation also has a critical point with a critical speed 
of 



c t 



1 + 3C 
-3-C" 



(21) 



Unlike equation (|17|), solutions cannot cro ss this critical 
point, since there dv/dz — > 00 (but see §3.2). 

We expect dp cx t/dz < 0, so solutions always accelerate 
(decelerate) for v > c-\ (v < Cf). Since c-f only exists for 
( < —1/3, solutions with ( > —1/3 always accelerate. In 
that case equation ( fl7| ) implies that for v > c+ (v < c*) the 
jet is expanding (contracting) in the r-direction. Since c* 
only exists for £ < 0, solutions with C > always expand 
sideways. 

If, on the other hand, ( < —1/3, two branches of so- 
lutions exist: (a) solutions which are injected with v > 
cj, which always accelerate and go through a nozzle at 
v = c* > c-f, and (b) solutions which are injected with 
v < Cj, which always decelerate. Thus, at sufficiently large 
dista nces from the core for gravity to be negligible (see 
§3^), highly anisotropic solutions with have to be 

injected at relativistic velocities to be accelerating, since 

— > 1 as £ — > — 1. This corresponds to the right branch 
of the dashed solutions pl otte d in Fig. [2] (which includes 
the effects of gravity, see §3.2). 

It is instructive to look at the case of pure anisotropic 
relativistic particle pressure again. We define the pressure 
anisotropy as ( p = (p|| - pi_)/{pi_ + py). If we fix C P by 
some rea rr ange ment process as we did for the magnetic 
field in § 2.1.1 (which might occur, for example, if there 
is coupling between magnetic field and turbulent pressure 
as suggested by B95), the behavior is very similar in the 
sense that accelerating solutions for C p > —1/3 (i.e., for 
p± > 2pii) have to be injected at super-critical velocity 

«>c tp =y(i + 3c P )/(5-c P ). 

If we simply let the components of p evolve adiabat- 
ically (without rearrangement), we arrive at a different 
critical velocity, c+ = \/(5 + 7( p )/(9 + 3£ P ) (note that 
C is not constant in this case), which exists only for 
Q p > —5/7. The solution once again has a nozzle at 
c* p =y(2 + 2C P )/(3 + Cp). Since c* p > c tp for £ > -5/7, 
and since dR/dz < for v < c± p , solutions injected with 
v > c-f must accelerate to satisfy pressure balance, which 
means that £ p increases with z, reducing c^ p and thus mak- 
ing the flow more super-critical (i.e., once above the critical 
point, the solution moves away from it). 

3.2. The Effects of Gravity on the Sonic Transition 

As seen in the previous section, a solution for £ > —1/3 
that starts out with v < c-f will always decelerate in the ab- 
sence of gravity. However, as is the case in the solar wind, 
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gravity can actually help a flow go through a critical point. 
We thus consider M > in this section. 

We can go through the same arguments as in §3.1. The 
critical speeds are still given by equations and (pl|), 
but now the critical conditions are different. At c*, equa- 
(14) gives 



tion 



[4p'„ (2 - 2 7ad ) -2(1 + C) tf'x*: 
+ (^ c 2 + V* + 2t^) ^ 



R~dz 







(22) 



instead of dR/dz = 0. Since £ > — 1 and 7 a d > 1, we 
can infer that dR/dz > at z*, i.e., there is no 'geomet- 
ric' nozzle at z+ anymore. The solution can always adjust 
dR/dz to satisfy this condition, thus the first critical point 
c* becomes irrelevant. 

Inclusion of the gravity term changes equation ( p0| ) to 



r 2 [(3 + 0« 2 + (i + 3C)] 



dv 
vdz 



(i + C) 



(3 + O^J 



(23) 



Now solutions can cross the critical point C|, since the 
right hand side of equation (E3) vanishes at 



3 + C\ 2r. 



(24) 



If t> ^ c+ at Zj, the solution must follow dv/dz = at 
that point. This is true for all £. Since that is the only 
zero of equation (|23|), we can therefore conclude that so- 
lutions accelerating at any z > z^ will accelerate for all 

Z > Zf. 

Solving the equations for dR/dz instead gives 



[(c 2 

=-(* 



1 



(3 + C) (i 

(C- 



0* 

' z 



1) 



•0] 

2rv 



(25) 



which also has a critical point at cj . For u = c+, the right 
hand side of this equation only vanishes at z+. In that 
case, dR/dz remains finite. For all other solutions (i.e., if 
z =/= z-f when v = C|), we must have singularities in both 
dv/dz and dR/dz. The singularity in dv/dz is evident 
from Fig. @ and from equation (p3|); pressure balance then 
requires that dR/dz must have a singularity of opposite 
sign, since dp ex t/dz is assumed to be finite, i.e. p C xt is 
continuous. 

We have numerically integrated equation ( p3| ) for two 
representative cases (£ = —0.9 and £ = 0, £ = 2) and plot- 
ted them in Fig. 0. Solutions are qualitatively different for 
C < -1/3 and ( > -1/3: 

• For £ < —1/3, there is one accelerating transonic 
solution, given by the condition in equation 
shown in the upper panel of Fig. |^ as a thick solid 
black curve. This is also the only solution acceler- 
ating for all z. As in the case of a regular adiabatic 
wind (Parker 1958), there also exists a decelerating 
transonic solution. Regions where solutions contract 
in the r-direction (i.e., where dR/dz < 0) are shown 



as hatched areas. There are four more branches of 
solutions. Two branches are double- valued (shown 
as dashed curves in Fig. ||). The left branch of those 
solutions can be rejected since those solutions only 
exist for z < Zj. For an accelerating solution on the 
right branch to exist, it must be injected with v > cj_. 
This corresponds to the solutions discussed in § jO| 
for which gravity can be neglected. 

The remaining two branches are solutions that are 
always sub- or supercritical (plotted as thin solid 
black curves in Fig. ||). The sub-critical solutions de- 
celerate for large z and always stay sub-relativistic. 
They are uninteresting as possible candidates for rel- 
ativistic jets. The supersonic solutions decelerate for 
z < z\ and accelerate for z > z\ . These solutions cor- 
respond to the super-critical solutions mentioned in 
§ |3.l[ . They always expand in the sideways direction. 
As we let £ — ► —1, z\ — > oo. This is not necessarily 
an indication that no solution is possible for -1, 
since for those cases the critical speed is very close 
to 1, thus the solution can attain large T. Further- 
more, as we saw above, the solution is expanding 
even before it goes through z+. We can thus have 
a regular (though sub-critical) accelerating jet even 
forC«-l. 

• For C > —1/3, there is only one branch of solutions, 
all of which start out decelerating, shown in the bot- 
tom panel of Fig. |2[ As the solutions reach z\ they 
begin to ac celer ate and behave the same way as de- 
scribed in §3.1. Since Z| moves inward for increas- 
ing C, this is no handicap. For ( > —1/3 we have 
Zj < 8r g /£ from equation (|24|), which, for reason- 
able values of £, is well in the regime where rela- 
tivistic corrections become important and inside the 
region where we expect the injection to occur. All 
of these solutions have positive sideways expansion 
dR/dz > for all z. 

The transonic solution for £ < —1/3 has some addition- 
ally nice features: Since we know z+, we can relate the 
jet cross section to the total jet power Lj at z\. Assum- 
ing that the jet is still magnetically dominated at z-\, the 
kinetic luminosity of the jet is 



L t = 7 ri? t 2 r t 2 w t 2l7 / _ L1 



(26) 



The external pressure at z-f is p e xt+ and so 



Ri 



\ 



LA 5 



i+C 
i-C 



(i + 



(27) 



While we do not know j? ext ± , for most parameter choices 
Lj is very nearly equal to L^, which we have a reason- 
ably good handle on from an observational point of view. 
Furthermore, we can estimate the jet width at observable 
distances and scale the solution back to Zf , which gives us 
an estimate of P] and thus U'±*. This in turn will allow 
us to determine the original matter loading of the jet from 
estimates of the terminal Lorentz factor T^. 
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2. — Inner solutions for £ = 2 in the vicinity of the critical 
I in the absence of dissipation but including gravity for the 



Fig 
point 

two possible cases, £ < —1/3 and £ > — 1/3. In the first case (panel 
a), shown for £ = —0.9, the transonic solutions are thick black lines, 
sub- and super-critical solutions are thin black lines. Dashed lines 
show double-valued solutions (of interest is only the upper right 
quadrant for 2 > Z* , v > c-j- ) . The critical values c* , Z* are shown as 
dash-dotted lines, c* is plotted as a grey dotted line. The condition 
dR/dz = is shown as a dash-tripplc-dotted line and regions where 
dR/dz < are shown as hatched areas. In the second case (panel b, 
for £ = 0) only super-critical solutions are possible. All solutions are 
expanding sideways and all possible initial values lead to acceptable 
solutions. We have once again indicated the location of z\. 



4. SOLUTIONS IN THE SELF-SIMILAR REGIME AND 
ASYMPTOTIC SOLUTIONS 

For an already relativistic jet in the 'self-similar' range 
Zf <C z -C ^inertia (where Zincrtia is the location where 
p'c 2 — 2U'±) a self-similar solution can be found. Equa- 
tions (§) and (jl|) give T oc R~^>/^ and with equation @ 
we have T oc R. 

Under the conditions of equation (|]), the pressure bal- 
ance condition gives 



T oc 7? OC Pext 



-1/4 



(28) 



the same as in the case of isotropic particle pressure con- 
sidered by BR74. For future reference we define the accel- 
eration efficiency 



rflnT 
dlnpext ' 



(29) 



thus for this simple case r\ = 1/4. 

If we ad opt the less limiting restriction \i\ / ^2 = Hz/ Hi 
(see § |2.1.l| ) instead of equation (||), we can find powerlaw 
solutions in three limiting cases: 

(a) For (5 = 0, the solution is given by T oc R^^ 2 /^ 1 oc 



Poxt 



M2/4/J1 



, which c an be very efficient for fj,i <C — //2 
(as pointed out in § |2.1.l| , one would generally expect 
that fi2 < < fix. 



(b) For the solution is given by T oc R 



-a»2/mi 



oc 



Poxt 

77 < 1/2 



1/(2^/^2-2) 



which has a limiting efficiency of 



(c) For (Ril the solution is approximately the same as 
case (a). 

If 5 > 0, the solution will in general approach solutions 
(b) or (c) (i.e., C ~~ > ±1)- For 6 < it is possible that the 
solution approaches a finite terminal Lorentz factor and 
zero opening angle if radial tension cancels the pressure 
due to and U^. Fig. || shows the different regimes. 
Note once again that these solutions exist only for r> 1 
and pc 2 < (/i. For all other cases the coefficients hi, p,2, 
/X3, and Hi are not constant. Note that we would generally 
expect that p\ ~ —p-2, since otherwise the re-arrangement 
process would be acting preferentially for changes in ge- 
ometry in one specific direction, which seems arbitrary. 
Thus, these results reduce to the well known r/ ~ 1/4. 

We can look for solutions in the presence of dissipation 
of magnetic energy. We now have to consider equation 
(|T|). We assume that the energy goes completely into 
relativistic particles, thus energy conservation implies 



dz 



dissipation 



I dU' 
3 dz 



(30) 



dissipation 



The particle energy can subsequently be radiated away as 
isotropic radiation. As long as p' <C U'±, the radiative 
case is no different from the non-radiative case, since the 
adiabatic term in the particle pressure does not contribute 
to the dynamics. 
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Fig. 3. — Regions in parameter space for the most general choice 
of parameters fii possible in the self-similar region. We have indi- 
cated the various limits of f at sufficiently low external pressure. 
Along with £oo we have indicated the limiting efficiency achievable 
when £ has reached is limiting value. For — /X2//J1 > b the efficiency 
rj can become very large. For <5 < 0, it is possible that the jet stalls 
due to radial tension. Note also that for 8 < there is a minimum 
value for £ at the injection below which no solution is possible, given 
by Cmin > (1 + S)/(S — 1). Even if £0 > Cmin; the solution will con- 
tract and decelerate if 6 is below a certain value. We have plotted 
this parameter range for a jet that starts out isotropically (i.e., with 
£ = —1/3) as a grey area. For smaller £, this area will become 
larger. 

A powerlaw type solution is once again possible only 
if the enthalpy is ne glig ible compared to the magnetic 
energy density (see §4.2). In that case the solution in 
the self-similar range is given by T oc Pext -7 ' with 77 = 
l/[4 + 6A(3 + 5C)/(3 - 3C 2 - 2A - 6A0] < 1, with rj = 1 
at C = — 1. We have plotted 77 as a function of £ and A in 
Fig. ||. For £ < —3/5 the efficiency is larger than in the 
case without dissipation, for £ > —3/5 it is smaller. The 
limiting efficiency that can be achieved in such a jet un- 
der the condition that dR/dz > is given by 77 < 1/2, as 
A — > 0. This happens because the dissipative process can 
convert energy in the parallel field component Uu (which 
does not enter eq. pljfl) into particle pressure, which must 
be taken into account in the energy balance. Also shown 
are areas in £-A-space where the conditions dT/dz > 0, 
dR/dz > are not satisfied. 

4.1. Opening Angles and Causal Contact 

Since a jet is generally defined as a collimated out- 
flow, we can ask under which conditions the solutions 
from above are actually collimated. The collimation 
condition dR/dz < 1 translates to a pressure gradient 
—d\np ext /d\nz — £ < 4 as long as the jet is magnetically 
dominated, the same as in the particle pressure dominated 
case. The presence of dissipation can alter this value. Gen- 
erally, the collimation is increased by dissipation, since the 
sideways pressure is reduced, thus the jet does not need to 
expand as much. In the coasting phase, where the jet is no 
longer accelerating, this condition changes to £ < 3+C < 4. 

Given the solutions from above, we can investigate 



the ratio of the opening angle a Q to the beaming angle 
o?b ~ r _1 . In the absence of dissipation we can write 



dR f/2-1 
ao/ab = I ~r °c z 1 
dz 



(31) 



independent of £. Thus, for steep pressure gradients £ > 2, 
the opening angle will grow faster than the beaming angle 
and will thus always become larger even if it starts out 
being smaller. For shallow pressure gradients £ < 2, the 
situation is reversed, i.e., the beaming angle will eventu- 
ally become larger than the opening angle. The presence 
of dissipation changes this behavior qualitatively: the ra- 
tio cko/ab now depends on both A and £, as illustrated in 
Fig. § 
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Fig. 4. — Contour plot of acceleration efficiency r; = 
— d In r/d In pext for radiative dissipation and different values of f 
and A. The grey areas indicate parameter values for which the jet 
is either not collimated (light grey with dark grey border) or deceler- 
ates (dark grey with light grey border). The hatched areas indicate 
regions where the ratio a /cfb of opening angle to beaming angle 
grows with z for given values of < £ < 4. 



This has consequences for the appearance of the jet, 
since the effective beaming angle is given by the larger of 
the two angles. Under the assumption that the jet is al- 
ways collimated, the opening angle in the coasting phase 
will always become smaller than the beaming angle, since 
the jet does not accelerate anymore. This could have im- 
portant consequences for the morphology of superluminal 
sources: if the beaming angle were smaller than the open- 
ing angle, one might expect to see larger jet misalignments, 
or lose the jet morphology altogether. The appearance 
would become sensitive to the emissivity and local Lorentz 
factor as a function of position across the jet cross section. 

Jets that expand too fast will eventually lose causal 
contact with their environment. This happens when the 
Alfven crossing time of the jet becomes larger than the 
expansion time (in the comoving frame), i.e., when 
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(where we approximated ^Aifven ~ c) or 



R > 



IT 



(33) 



This corresponds (up to the factor £) to the criterion when 
a a > ah- Thus, for £ > 2 (in the absence of dissipation) 
the jet will eventually lose causal contact with its sur- 
roundings (see Fig. [| for values of £ and £ where this is 
the case). As mentioned in §|], a quasi ID treatment is no 
longer possible, since the internal pressure balance is now 
regulated by shocks traveling inward from the jet walls. 
After the jet reaches the terminal phase, it will re-gain 
causal contact, since the opening angle will continually 
decrease (assuming the jet is still collimated). 

4.2. Equipartition 

Constant pumping of magnetic energy into particle pres- 
sure can lead to equipartition between U' and p. We can 
use the self-similar solution to estimate r cq , where the ac- 
cumulated particle pressure surpasses the magnetic energy 
density (including effects of adiabatic cooling on the accu- 
mulated particle pressure, where we assume that it behaves 
as a relativistic gas, which gives an upper limit on the re- 
sulting pressure). Thus, for J> r eq the solution might 
be altered. For some parameter values p never reaches the 
level of U'±. In that case we estimate the asymptotic ratio 
(p' /U' ±) qo . Figure|| shows the results of those estimates. 
For large £, equipartition can be reached quickly, thus, 
unless the energy going into particles is subsequently ra- 
diated away, the assumption that the particle pressure be 
negligible compared to the magnetic field energy density 
will be violated beyond r oq . 
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Fig. 5. — Value of T for which the pressure accumulated by non- 
radiative dissipation reaches equipartition with U' ±. Neglecting p 
is no longer justified beyond that T. For £ < —0.6 the pressure 
never catches up with U'±, in that case we plotted contours of the 
limiting ratio p'/U' ± (dashed line^). The hatched and grey regions 
correspond to the regions in Fig. 

We define the energy distribution function as 

r72 



where 7 is the Lorentz factor according to a particle's 
random motion, measured in the comoving frame, and 
7i *C 72 are the lower and upper spectral cutoffs. Since we 
assume that the magnetic field is dominating the internal 
energy budget, synchrotron losses can be very strong, pro- 
vided the particle energy spectrum is flat enough so that 
most of the energy is at high particle energies (i.e., s < 2). 
In that case we can expect most of the energy to be ra- 
diated away immediately and the corresponding electrons 
will lose most of the inertia, thus the dissipated energy 
will not lead to a build-up of particle pressure. If, how- 
ever, synchrotron losses are weak compared to dissipation 
(e.g., if the spectrum is too steep, or if synchrotron self- 
absorption traps most of the radiation to inhibit cooling), 
the effects of particle pressure can become important, as 
demonstrated in Fig. ||. For a discussion of the observa- 
tional effects of the different radiative scenarios see 



5.1 



/( 7 ) = F o7 - s , / /(7)d 7 = ri, 



(34) 



4.3. Full Solutions 

We can solve the full equation ( |l4| ) in the regime T ^> 1 , 
i.e., for relativistic jets. As mentioned before, the pressure 
balance condition leads to a simple algebraic equation in T 
and R. In the absence of dissipation and gravity, equation 
(H3) is in fact another algebraic equation relating T and i?, 
thus, the two equations can be solved for r(p cxt ) using a 
numerical root finder. Apart from reproducing the scaling 
behaviors established in §0, this will enable us to determine 
the terminal Lorentz factors and the length scales over 
which the transitions between different dynamical phases 
occur. Furthermore, we can use the full dynamical model 
to investigate the evolution of such observational quanti- 
ties as polarization and synchrotron brightness. 

In the absence of dissipation, the terminal Lorentz fac- 
tor that can be reached with such a jet is simply 



lim F^F f ^ 2 + 2 ^ + 4 ^ ), (35) 



where subscripts denote quantities evaluated at some 
arbitrary upstream point. 

This simple solution is no longer possible in the presence 
of dissipation, which introduces a sink term into equation 
(|l3|). As a result, we have to use equation (|^) instead. 
Once the energy has been converted into particle pres- 
sure, it can be radiated away as isotropic radiation, which 
will not affect the dynamics of the jet any further (as- 
suming that p is dynamically unimportant). If the energy 
is stored as particle pressure, the pressu re c ould eventu- 
ally become dynamically important (see § |4.2| ). Until that 
happens, though, the two solutions are identi cal . The ter- 
minal Lorentz factor is always reduced (see §^|), but the 
acceleration efficiency can be increased for £ < —3/5 (see 
§^). We have plotted the solution for the radiative case 
(the one case solvable analytically) in Fig. ||. 
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— pure magnetic jet (self similar case) 
. ■ ■ cold matter (2U perp =40p c 2 ) 

— cold matter, dissipation for A=0.05 
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Fig. 6. — Analytic solutions in the limit v ~ 1 for f = 0: a) self 
similar limit, valid for very large U± 2> p'c 2 ; b) solution with cold 
matter, Uj_ /poc 2 = 20, i.e., Too = 41; c) radiative dissipation with 
A = 0.05, otherwise same parameters as b); d) same parameters as 
c) but including radiation drag. Note that the dissipation of en- 
ergy alters both the acceleration efficiency and the terminal Lorentz 
factor Too. 



4.4. Radiation Drag 

The presence of ultra-high-energy particles in AGN jets 
suggests that inverse Compton (IC) enhanced radiation 
drag might be dynamically important. While O'Dell 
(1981) initially suggested that pair jets might be acceler- 
ated by the Compton rocket effect, Phinney (1982) showed 
that it is hard to accelerate a plasma beyond fairly modest 
Lorentz factors by radiation pressure without a continuous 
source of particle acceleration to offset the strong IC cool- 
ing of the plasma. Furthermore, if the radiation source 
is not point-like, the terminal Lorentz factor is limited 
by the solid angle f2 the radiation source subtends. On 
the other hand, radiation drag can hamper the bulk accel- 
eration of plasmas containing relativistic particles in the 
presence of a radiation field, if those particles are continu- 
ously reheated to overcome the IC losses. The dissipation 
mechanism discussed above could provide such reheating. 
Here we will consider the effect of radiation drag in the 
simplest possible prescription. 

We assume that the IC cooling and the dissipational 
heating time scales are short compared to the adiabatic 
timcscalc. If this is not satisfied, the influence of radiation 
drag will be reduced. We can then expect dissipational 
heating to nearly balance IC losses in a near equilibrium 
situation. Thus the amount of IC drag is controlled by 
how much dissipation there is. For this approximation to 
be valid, IC losses must dominate the loss processes of the 
particles, i.e., the radiation energy density C/' ra d in the 
comoving frame must be large compared to the magnetic 
field energy density U' (for large enough T this is always 
going to be the case, since the external field will be Doppler 
boosted). Finally, we assume that (/? 2 7 2 ) 3> 1, where 7 
is the particle Lorentz factor in the comoving frame. This 



sets an upper limit of 

f/'rad < 6 x 10 3 ergs cm" 3 
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on the comoving radiation energy density (otherwise IC 
cooling would have lowered the upper spectral cutoff to 
72 ~ 1). These assumptions allow us to eliminate [/' ra d 
from the equations, since the drag term and the cooling 
term are both proportional to t/' ra d- In a sense, then, we 
are presenting an upper limit on the importance of IC ra- 
diation drag over large length scales. It has to be kept 
in mind, though, that drag can be much more important 
in non-stationary situations (like, for example, in shocks), 
which are beyond the scope of this paper. 

We assume the jet is moving through a radiation field 
that is locally isotropic in the lab frame. Following Phin- 
ney 's treatment (1982), we can calculate the loss rate and 
the force density due to IC scattering in the comoving 
frame and then transform back to lab frame to find the 
additional term for equation (|l3|). We find that radiation 
drag always decreases both the acceleration efficiency 77 and 
the terminal Lorentz factor by moderate amounts. It 
does not, however, introduce qualitatively new features. 
To demonstrate this, we have plotted a solution including 
radiation drag for otherwise identical parameters in Fig. 

5. DISCUSSION 

5.1. Tradeoff Between Dissipation and Acceleration and 
Synchrotron Brightness 

In the following section we will investigate the observa- 
tional effects of the jets we have introduced in this paper. 
A highly dissipative jet will radiate away a large fraction 
of its internal energy along the way before reaching the 
terminal Lorentz factor Too, while a non-dissipative jet 
will convert all its internal energy into kinetic energy flux. 
Since the jet will ultimately terminate and reconvert its 
kinetic energy flux into random particle energy when it 
slams into the surrounding medium, the ratio of kinetic 
luminosity (which could be estimated based on the energy 
input into the lobes, based on the source size and its age) 
to the radiative luminosity L^i SS (z) (i.e., the integrated 
luminosity of the jet before reaching the terminal shock) 
should give us some indication of the importance of dissi- 
pation. 

We have already seen in § [4.3| that the presence of dissi- 
pation can lower , thus lowering the kinetic energy flux 
at the terminal shock, Loo (dominated by cold kinetic en- 
ergy flux), and the produced hot-spot luminosity. A given 
fraction b of the terminal luminosity bL^ will be radiated 
away, which can be estimated from the hot-spot and co- 
coon luminosity (calculating b is, of course, a highly non- 
trivial matter), giving us a handle on A. We have plotted 
the ratio 



(37) 



in Fig. 0. To make that plot, we chose parameters such 
that in the absence of dissipation the jet would reach 
T 00 (A = 0) = 2U' u /(p' Q c 2 ) = 100. If we had chosen 
a larger (smaller) value of this parameter, the lines in the 
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plot would move down (up), since depends non-linearly 
on both A and 2U' /(p' c 2 ), so this plot is only a rep- 
resentative one of a family of similar plots for different 

r 00 (A = o). 




Fig. 7. — Plot of the ratio Aiiss/^oo, indicating how much of 
the energy carried in the jet is radiated away and how much will 
reach the terminal shock (this energy can be used to heat the hot 
spot and to inflate the radio lobes). This plot was constructed for 
Tooo = Too (A = 0) = 100. Lower values of will shift the lines 
in this plot upwards. Note that this plot was constructed using the 
assumption that all the dissipated energy is radiated away on the 
spot. 

Another question is what the spatial distribution of the 
jet emission is, since there are several competing effects: 
the optical depth to self absorption (and inverse Comp- 
ton up-scattering), Doppler boosting and opening angle 
(the larger of which will determine the opening angle of 
the cone into which most of the radiation goes), and of 
course the dissipative power of the jet itself [depending on 
A and T(z), R(z)]. This question can be asked with re- 
spect to the frequency integrated brightness I or the spec- 
tral brightness l v . If we assume that all the dissipated 
energy is radiated away on the spot, we can calculate the 
local dissipation rate, which must then be equal to the lo- 
cal frequency integrated jet emissivity f in the comoving 
frame. This assumption depends on the injected particle 
energy spectrum. If the spectrum is flatter than s = 2, 
most of the energy is in the high energy particles. In that 
case, synchrotron cooling can be efficient enough for our 
assumption of on the spot radiation to be effective. If, on 
the other hand, s > 2, most of the energy is in the low en- 
ergy particles, synchrotron radiation will not be efficient 
(unless 71 is very high, in which case the injected spectrum 
would rapidly cool to a quasi mono-energetic distribution), 
and the jet will accumulate particle pressure or radiate by 
other means (note that Compton cooling would be equally 
insufficient to balance heating in this case). 

To investigate the first case we will set s < 2. Given 
a viewing angle we can then determine the observed total 
intensity J, given by 



where 8 is the angle between line of sight and jet axis and 
D = [r (1 — v cos $)] 1 is the Doppler factor. This expres- 
sion takes relativistic beaming and the relativistic correc- 
tions to foreshortening into account. One might expect 
that the integrated brightness peaks at a certain distance 
from the core, since D is strongly peaked at T ~ 1/8. How- 
ever, in our prescription the dissipation drops off too fast 
for this effect to be important. The main difference in the 
brightness evolution is that a dissipative jet has a different 
efficiency r\ and generally expands less rapidly in the side- 
ways direction (and will reach a smaller terminal Lorentz 
factor Too). This effect will become important once the jet 
has reached and only for A large enough to significantly 
alter the dynamics (A > 0.1). We have plotted / as a func- 
tion of z for different values of A, £ = 0, U'o = 20 p' c 2 , 
and a viewing angle of 9 = 10° in Fig. ^|, arbitrarily nor- 
malized to I(A = 0.01) to increase dynamic range (the 
brightness decreases by many orders of magnitude along 
the jet). As is obvious from the plot, for small A, only the 
overall normalization of / varies with A, whereas for large 
enough A, the brightness distribution itself changes shape 
due to the altered dynamics. Also shown are the Doppler 
factors D for the different parameters, which are primarily 
responsible for the different shapes. 

The situation can become more difficult in the opposite 
case, i.e., if most of the radiation is trapped (e.g., by syn- 
chrotron self-absorption, which implies a steep spectrum, 
s > 3) or if the deposited energy is simply not efficiently 
radiated (e.g., if 2 < s < 3). In that case the emission 
of dissipated energy could be delayed, leading to a rela- 
tive brightness peak downstream. To briefly investigate 
this possibility, we assume the latter case, i.e., 2 < s < 3, 
which is not an unreasonable choice for AGNs (see § |5.3| , for 
example). We assume that the energy flux vF v peaks at 
high energies v v : either at the spectral break of Aa ~ 1/2 
expected in a scenario in which high energy particles are 
constantly re-injected, where the spectral index a is given 

by 
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or at the spectral cutoff produced by synchrotron and 
IC cooling (in the absence of a strong break). The po- 
sition of i/p depends on adiabatic effects, radiative cool- 
ing, and heating due to dissipation. The spectrum will 
be self-absorbed at low frequencies, which generally leads 
to an observed spectral index of a ~ —5/2 (for an ex- 
act treatment of the spectral shape at the self absorption 
turnover see De Kool, Begelman, & Sikora 1989). If we 
take 2<s<3orl/2<a<l, the self absorbed part con- 
tributes a negligible fraction to the total brightness and 
most of the energy is emitted at the high end of the spec- 
trum. 

In this case it is impossible to calculate the brightness 
analytically as a function only of T, R, and z. Rather, 
one can numerically integrate the evolution equation of 
the peak frequency under adiabatic cooling, dissipative 
heating [we assume a self-similar transfer of energy from 
magnetic field to the particles such that (d'y / dz)\ diss cx 
7 (dlf /dz)\ diss ], and synchrotron cooling on the basis of 
the solution T(z) given above. We can estimate the run of 
I by scaling it with the brightness at v v , taking account of 
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relativistic beaming and aberration. If we choose to nor- 
malize the intensity curves as we did in the previous case, 
we can get around fixing the absolute normalization of U', 
since it only enters linearly into the intensity and will thus 
cancel out upon normalization. We used U'± = 20 p' c 2 , 
along with s = 5/2 and plotted the frequency integrated 
intensity in Fig. with otherwise the same parameter val- 
ues as in Fig. g, once again normalized to / for A = 0.01. 
Note that for large values of A our assumption that the jet 
be magnetically dominated and that particle pressure be 
negligible can break down (see Fig. ||) , so curves with high 
A are to be taken with a grain of salt. Nevertheless, it is in- 
teresting to note that the intensity drops less rapidly with 
z for larger values of A, corresponding to the delayed emis- 
sion mentioned above. The shapes of these curves depend 
only weakly on s. The different bends in the curves stem 
from the evolution of the Doppler factor (shown in Fig. [|), 
and the evolution of the peak frequency. The slopes of the 
curves are produced mostly by the evolution of the lower 
cutoff frequency 71 (see eq. p4| ) , which enters the expres- 
sion for the synchrotron brightness through the powerlaw 
normalization of the particle distribution, and by the evo- 
lution of Lorentz factor and jet radius (entering through 
the particle density and the integration of the emissivity 
across the jet). 
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Fig. 8. — Plot of the frequency integrated brightness for a ra- 
diative jet (i.e., all the dissipated energy is radiated away on the 
spot) as a function of 2 for different values of A, £ = 0, f = 2, 
U' Xq/p'o = 20, and 6 = 10°, arbitrarily normalized to the intensity 
curve for A = 0.01 to increase the contrast (thick black curves). Also 
shown are the corresponding Doppler factors D as thin grey curves. 
The shapes of the individual intensity curves are mainly determined 
by the variation of D. Note that a significant observable effect is 
achievable only for A > 0.1. 



10 4 



10 j 



§ 100 



< 

< 



10 



A=.004 

A=. 010 

. . . A=.027 
.....A=.074 
_,.._A=200 
__ A=.541 



/ 

/ 

/ 

/ 



* 



10 



100 



10 5 



10 4 



10 5 



z/z Q 



Fig. 9. — Plot of the frequency integrated brightness for a 
marginally non-radiative jet (only a dynamically small amount of 
the dissipated energy is radiated away on the spot) as a function 
of z for different values of A, ( = 0, U'o/p' c 2 = 20, £ = 2, and 
= 10°, arbitrarily normalized to the intensity curve for A = 0.01 
to increase the contrast (the absolute intensity drops by many or- 
ders of magnitude). Note that these curves do not depend on the 
absolute value of the magnetic field, which cancels out due to the 
normalization. 



5.2. Polarization 

While the degree of polarization is highest for homo- 
geneous magnetic fields, jets with tangled or disorga- 
nized field can exhibit a net polarization if there is a net 
anisotropy in the field (Laing 1980). The measured polar- 
ization will depend not only on £, but also on the viewing 
angle 8 and the bulk Lorentz factor T. The polarization of 
radiation from a powerlaw distribution of electrons with 
index s in a region of homogeneous field is given by 



n : 



I±-Ii 



5 + 1 

3 + 7/3 



(40) 



where / is the intensity at a given optically thin wave- 
length. To calculate the integrated polarization, we aver- 
age across the jet. To do this we decompose the radiation 
into polarization along the jet axis and perpendicular to 
it. Furthermore, we assume that field directions are dis- 
tributed among all solid angles and introduce a weighting 
function that distributes the field orientations to the re- 
quired anisotropy, 



w(d) = sintf K 



(41) 



where 1} is the angle between the field and the z-axis and k 
is determined by the anisotropy. We can solve for k under 
the condition that (B«) = (1 + C)/(l - ()(B± 2 ): 



1 + 3C 
1 + C 



(42) 



We correct the viewing angle for relativistic aberration, 
which has a significant impact on the observed polariza- 
tion, since the average polarization will go to zero for a jet 
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seen head on. Furthermore, the angle between line of sight 
and magnetic field is important in determining the rela- 
tive brightness of a region. Figure [l(] shows the predicted 
polarization for the cases shown in Fig. |(| a spectral in- 
dex of a = 0.5, and a viewing angle of 9 — 10°. Since 
the anisotropy of the jet is fixed, the variation in the po- 
larization 11(0) is solely caused by changes in the viewing 
angle due to relativistic aberration. Generally, the polar- 
ization will be perpendicular to the jet axis if ( < —1/3 
and parallel if £ > —1/3. As long as equation @ holds, an 
extremum in 11(9) will be present and it should indicate 
the position where T = l/sin0, i.e., where the viewing 
angle corrected for aberration is 9' = 90°. Note that this 
polarization is averaged across the jet. In order to com- 
pare these predictions to actual measurements a relatively 
small correction for the emission weighted averaging across 
the jet at different angles must be made. The qualitative 
predictions of this section should be unaffected by that 
caveat. 



pure magnetic jet (self similar case) I 

cold matter (2U perp =40p c 2 ) 
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Fig. 10. — Optically thin polarization It as a function of z for the 
same parameter values as in Fig. M and a viewing angle of 9 = 10° . 
We chose a spectral index of a = 0.5. The variation in II is solely 
due to changes in the aberrated viewing angle along the jet. 

5.3. Applications 

Finally, we sketch out some applications for this model. 
The obvious candidates for jet models are radio galaxies 
and all related jet powered AGNs. The best studied ex- 
ample is M87, since it is the closest unobscured source 
(though by no means a particularly powerful one). The 
central black hole has a relatively well determined mass 
of M ~ 2.4 x 10 9 M Q or r g ~ 3.5 x 10 14 cm. The ki- 
netic power of M87 is estimated to be of order Lkm ~ 
10 43-44 ei . gssec -i (Reynolds et al. 1996, Bicknell & Begel- 
man 1996). Taking the new proper motion measurements 
based on HST data (Biretta, Sparks, & 1999) at face value, 
the terminal bulk Lorentz factor likely falls into the range 
6 < < 10 with a viewing angle of 9 ~ 20°. [Note that 
VLA measurement show slower proper motions (Biretta, 
Zhou, & Owen 1995) and the relationship between the ob- 
served pattern speeds and T is not known. There are, 
however, no direct measurements of bulk motions, so for 



lack of better knowledge we will use the larger values ob- 
tained from the optical data.] The average polarization 
at kpc distances (where the jet has likely reached termi- 
nal velocity) is roughly n ~ 10% parallel to the jet axis 
, which corresponds to £ ~ 0.3 with the numbers given 
above. 

Reynolds et al. (1996) showed that the jet probably con- 
sists of pairs rather than ionized gas, so the jet might ac- 
tually still be accelerating (if indeed a large fraction of the 
particles has 7 > 10. Note, however, that the presence of 
shocks, clearly visible as knots in all wavelengths, calls for 
a more sophisticated model which takes time-dependence 
and MHD instability effects into account. For a discus- 
sion of the nature of the shocks observed in M87 see Bick- 
nell & Begelman 1996.) Furthermore, there is evidence 
that the jet is not magnetically dominated on large scales 
[z > 100 pc) and that the magnetic field is actually some- 
what below equipartition (Heinz & Begelman 1997), which 
means in this context that the magnetic field must have 
been dissipated non-radiatively, accounting for the par- 
ticle pressure at larger distances. The observed spectral 
index at optically thin wavelengths is a ~ 0.6, which is 
at least in the correct regime for synchrotron radiation to 
be in efficient at radiating away the dissipated energy (see 

At a viewing angle of 20° the jet radius at knot A 
(z ~ 3 kpc) is R ~ 35 pc with an approximate opening an- 
gle of a ~ 0°.7, much smaller than the beaming angle of 
a b ^6°, thus the jet is narrow at least at VLA resolution. 
Note that entrainment of ambient material may be impor- 
tant before the jet reaches knot A (Bicknell & Begelman 
1996). Using knot D as a reference point does not change 
this analysis significantly, however. There are no good es- 
timates of the pressure gradient in the innermost regions of 
the M87 X-ray atmosphere (this situation should change, 
though, with the launch of Chandra). The pressure at 
large distances is approximately pism ~ 10~ 10 dyncm~ 2 , 
but Bicknell and Begelman (1996) argue that the pres- 
sure in the radio lobes is significantly higher than this, 
Pcxt ~ 1.5 x 10~ 9 dyncm~ 2 . For lack of better information 
we will assume the latter value and a smooth pressure gra- 
dient with £ = const. The radiative luminosity of the jet 
itself is L ra d > 3 x 10 42 ergssec -1 (Reynolds et al. 1996), 
which argues for A > 4 x 10~ 3 . If (as we speculated above) 
much of the particle pressure at large distances was indeed 
produced by dissipation, A could be much larger. 

The jet probably reaches terminal Too ~ 10 somewhere 
between VLBI scales and VLA scales, 2 x 10 3 r g ~ 0.2 pc < 
Zoo < 200 pc ~ 2 x 10 6 r g . For A small enough to be 
dynamically unimportant we can assume that rj ~ 1/4 
in the self similar regime. Arbitrarily setting zo ~ 10 r g 
gives 0.8 < £ <; 1.74. The central pressure po is then 
2 x 10 _4 dyncm~ 2 < po < 200dyncm~ 2 , which requires 
an rms magnetic field of 0.09 G < B'o < 90 G for pressure 
balance (using ( ~ 0.3). The initial jet width strongly de- 
pends on £: 8 r g < Rq < 9000 r g for an assumed zo = 10r g . 
The latter value is unrealistic (and inconsistent with the 
limits put on the jet width by VLBI observations, Junor & 
Biretta 1995), since most of the energy output of the disk 
into the jet will be provided close to the black hole. It is 
therefore most reasonable to assume that z^ ~ 2 x 10 3 r g 
and £ ~ 1.7. The total jet power implied by the num- 
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bers given is of order L ~ 2 x 10 44 ergs sec^ 1 , consistent 
with the estimate by Bicknell & Begelman (1996). Over- 
all it seems that this model is consistent with the observed 
properties of M87 to first order if we adopt a pressure gra- 
dient following £ ~ 1.7. 

The total isotropic energy output of E > 10 54 ergs of 
GRB 990123 (Bloom et al. 1999, Bloom et al. 1999) ar- 
gues strongly in favor of non-isotropic gamma-ray burst 
(GRB) scenarios, so jet models explaining the apparently 
beamed nature of these sources (e.g., Meszaros & Rees 
1997, Sari, Piran, & Halpern 1999) enjoy newly enhanced 
popularity. One of the standard scenarios for the energy 
sources of GRBs is a massive accretion event (either a 
neutron star - neutron star merger or the accretion of a 
neutron star by a black hole: Narayan, Paczyhski, & Pi- 
ran 1992), which leads to the formation of a disk after 
tidal disruption of one of the objects. Since neutron stars 
already display large magnetic fields, one might expect 
strong shear amplification of this field in the disk, leading 
to a scenario similar to what we described in Similarly, 
the hypernova approach (Paczyhski 1998, Woosley 1993) 
can produce a precollimated outflow, which might then 
evolve into a jet. Application of our model to GRBs is, 
however, not as straightforward as in the case of mature 
radio galaxies. This is because GRBs are highly time de- 
pendent (corresponding to the adolescent stages of radio 
galaxies): the lifetime of the jet (roughly of the order of 
the light travel time of the material) is shorter than the 
sound crossing time of the bubble the jet blows into the 
environment, so a pressure balanced solution as described 
above (which will be set up after the jet and the ambient 
material have equilibrated, i.e., after the jet has existed for 
a few sound crossing times) might not be a good approxi- 
mation. The investigation of jet dynamics in the context of 
GRBs and in the framework of tangled fields as introduced 
above will therefore be the subject of another paper. Here 
we simply wish to point out the benefits of jet models in 
general and an approach based on our model in particular: 

• Acceleration of GRB outflows by Poynting flux has 
the advantage that collimation can be provided not 
only by the external medium, but also by the field 
geometry itself (note that for 5 — the perpendic- 
ular component of the field does not contribute to 
the sideways pressure, so for £ ~ — 1 the jet can 
have a large Poynting flux yet orders of magnitude 
smaller sideways pressure than a particle dominated 
jet would have for the same energy flux). 

• One would not need to invoke shocks to produce 
emission in this context, since the internal dissi- 
pation of magnetic energy would provide a natural 
source of high energy particles and photons to pro- 
duce gamma rays (see Thompson 1994 for an ex- 
ample of how internal dissipation can power GRBs). 
The short term variability seen in GRBs could then 
be explained by inhomogeneities (e.g., variations in 
the inhomogeneity of the field) imprinted on the out- 



flow by variability in the central engine itself, which 
is expected to have time scales of the same order as 
the ones observed. 

6. CONCLUSIONS 

We have presented simple analytic models of jet acceler- 
ation. The jets are accelerated by tangled magnetic fields, 
with collimation being provided by pressure from an ex- 
ternal medium. This is a new approach, based on Begel- 
man (1995), as previous models of MHD jet acceleration 
concentrated mostly on large scale organized fields. Our 
analytical quasi-lD approach is limited to narrow jets, for 
which the opening angle is smaller than the beaming an- 
gle, coincident with the condition that the jet be in causal 
contact with its environment. We introduce an ad-hoc pro- 
cess that redistributes energy between perpendicular and 
parallel field to facilitate efficient conversion of internal en- 
ergy to kinetic energy. Without such a process, stationary 
jet acceleration by disorganized fields is impossible. In the 
absence of dissipation, the rate of acceleration achieved 
under this scenario is the same as in the case of relativis- 
tic particle pressure, T oc p xt -1 ^ 4 - We also find analytic 
solutions beyond the self-similar region, which enable us 
to calculate the terminal bulk Lorentz factors of such jets. 
In order for these jets to reach the observed T > 10, they 
initially must be magnetically dominated. 

We estimate the impact of dissipation of magnetic en- 
ergy on the dynamics of the flow by considering a simple, 
phenomenological prescription of the loss process. The 
presence of dissipation lowers the terminal Lorentz fac- 
tor Too and generally changes the rate at which the jet is 
accelerated (the latter effect is noticeable only if the dissi- 
pation rate is comparable to the adiabatic expansion rate). 
We also include the effects of radiation drag in the simplest 
scenario, which always lowers the efficiency r\ and Too. The 
amount of radiation drag in our model is controlled by the 
amount of dissipation replenishing the high energy particle 
pool but seems to be dynamically unimportant. 

We calculate the frequency integrated surface brightness 
for the extreme cases where all or very little of the dissi- 
pated energy is radiated away on the spot and find that, 
while the brightness drops off very rapidly in all consid- 
ered cases due to the expansion of the jet, values of the 
dissipation efficiency A ;> 0.05 can have a significant im- 
pact on the intensity as a function of z. In the marginally 
non-radiative case the buildup of particle energy leads to 
a slower decline in intensity with z for larger A. Finally, 
we apply this model to the prototypical radio galaxy M87 
and find that it is consistent with the observed properties. 
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APPENDIX 



A MORE REALISTIC DISSIPATION LAW 



The dissipation law we assumed in §2.1.2 was only one of many plausible ad hoc models. Since (for finite conductivity, 
i.e., beyond the limit of perfect MHD) reconnection will occur whenever there is field reversal on sufficiently small scales 
(which would certainly be the case in a highly tangled geometry), we would expect dissipation to occur even if there 
were no change in field geometry due to expansion or acceleration. In that case we might expect that the dissipation 
timescale is proportional to the time it takes a disturbance to travel a given characteristic length (e.g., the jet width) in 
the comoving frame, i.e., 



dU'i 



dz 



A , ^Alfvcn 

1 TR 



(Al) 



Once 



where the parameter A absorbs the effects of resistivity and all the unknown physics of the reconnection process 
again, it is straightforward to generalize to the case of different Aj for different components of the field. 

It is possible to solve the set of equations in the magnetically dominated case (i.e., setting p' + Ap' = 0), under the 
assumption that £ = const, and in the relativistic limit, r> 1. In this limit i>Aifvcn = 1, which simplifies the treatment 
significantly. We assume that the dissipated energy is radiated away immediately, which leads to a modified equation (14) 



A f 2 + 6 C 
TR V3- 3C 



„*\ dT . dR 
^ + 2 ^fdz^ 2 + ^Rdz 



0. 



The pressure balance equation (|T^) is also modified: 



(C-i) 



dr 



(3 + C) 



dR 



A 

TR 



Tdz Rdz 

where we used U' oc p cx t and equation (|l9|). This set of equations can be solved to give 

(3+5C)/(4+4C) 



\z) = T (z/z f 4 [l-A [(z/zo) 1 -^ 2 - l] } 



with 



A = 



A 



1 



r i? (l-£/2)(l + 3C) 



(A2) 



(A3) 



(A4) 



(A5) 



Subscripts denote quantities evaluated at some arbitrary upstream point zq. In the limit of A — > this solution appropri- 
ately reduces to the result without dissipation, i.e., T oc z^/ 4 . The solution can essentially take on three different behaviors: 
If £ > 2, the solution will asymptotically approach T oc z^ 4 , i.e., n = 1/4. If £ < 2, on the other hand, two different 
scenarios can occur: the flow can either approach a self-similar behavior with 77 =/= 1/4 (but still constant) in the limit of 
z > z , or, if A > (i.e., C > -1/3), the flow can actually stall, i.e., T -> for z -> z [(1 - £/2) (1 + 3C) /A - l] 1 /^"^ 2 ). 
In this limit, the magnetic field dissipates away too quickly to satisfy pressure balance and the jet must contract and 
decelerate to increase its internal pressure, thereby increasing its dissipation rate, which leads to a run-away process. 
Eventually, the massless aproximation will break down, in which case the Alfven velocity will drop, lowering the dis- 
sipation rate (furthermore, the particle pressure will gain in importance, eventually stabilizing the jet against external 
pressure, in which case the jet would behave as described by BR74). This process would produce an observable hot-spot 
(and possibly a shock) at a fixed distance. 



